Minimal model for active nematics: quasi-long-range order and giant fluctuations 
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We propose a minimal microscopic model for active nematic particles similar in spirit to the 
Vicsek model for self-propelled polar particles. In two dimensions, we show that this model exhibits 
a Kosterlitz-Thouless-like transition to quasi-long-range orientational order and that in this non- 
equilibrium context, the ordered phase is characterized by giant density fluctuations, in agreement 
with the predictions of Ramaswamy et al. [Europhys. Lett. 62, 196 (2003)]. 

PACS numbers: 05.70.Ln,87.18.Ed,45.70.-n 



Over the last decade or so, physicists have been 
looking for common, possibly universal, features of the 
collective motion of animals, bacteria, cells, molecu- 
lar motors, as well as driven granular objects P, 0. 
Among the emergent properties of these groups of 
"active" or self-propelled particles (SPP), distinctively 
out-of-equilibrium features have been found, such as 
the existence of long-range orientational order in two- 
dimensional "ferromagnetic" flocks of polar SPP 0, Q . 

Another set of striking intrinsically non-equilibrium 
properties have recently been predicted by Ramaswamy 
and co-workers 0, 0, They considered in particu- 
lar the case of apolar but oriented SPP and argued that 
such "active nematics" should differ dramatically from 
the usual (equilibrium) case 0. In particular, their ap- 
proach, based on the analysis of hydrodynamic equations 
derived from symmetry arguments, predicts that giant 
density fluctuations arise in the ordered phase of such me- 
dia. In , it is also hinted at the possibility of true long- 
range order and of a different isotropic-nematic transition 
out-of-equilibrium, but no definitive statement is offered. 
Resolving these issues is nevertheless crucial, especially 
in view of the predicted giant density fluctuations, and 
all the more so since, in polar SPP, the transition to true 
long-range order was shown to be discontinuous Q . 

In spite of the above-mentioned current surge of activ- 
ity in non-equilibrium systems, the giant density fluctua- 
tions predicted by Ramaswamy and co-workers have not 
been observed so far, and the nature of the nematically- 
ordered phase and of the transition leading to it have not 
been elucidated. Experimentally, relevant systems such 
as colonies of elongated cells |9t and ensembles of rod-like 
objects driven by vibration [lOlllllll^ have been studied, 
but with other issues at stake. On the theoretical side, no 
microscopic model has been proposed this Letter, 

we fill this gap, confirm for the first time the predictions 
of Ramaswamy et al, and investigate the nature of the 
isotropic-nematic transition in driven systems. We in- 
troduce a minimal model for active nematics, and show 
numerically that its isotropic-nematic transition in two 
space dimensions does not differ significantly from the 



equilibrium case: only quasi-long-range (QLRO) order is 
attained, with scaling laws compatible with those of the 
Kosterlitz-Thouless (KT) transition. Nevertheless, giant 
density fluctuations are clearly observed in the ordered 
phase: the standard deviation An of n, the average num- 
ber of particles in a given sub-system, is proportional to 
n and not to ^/n as expected in equilibrium. 

Our model is similar in spirit to the Vicsek model for 
polar SPP 0. In a typical driven-overdamped dynamics, 
identical pointwise particles move synchronously at dis- 
crete timesteps At by a fixed distance voAt. In two space 
dimensions and for uniaxial nematics — the case to which 
we restrict ourselves in the following — , each particle j is 
endowed with an orientation 6j and moves along 9j or 
9j + TT with equal probabilities. At every timestep, 9*~^^ 
is given by 0(Q*), the direction of the first eigenvector 
of the local tensorial traceless order parameter 
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where the average is taken over all particles k within the 
interaction range ro — 1 > vqAI, including particle j (in 
this paper, we use vqAI — 0.3). As for the Vicsek model, 
disorder arises from the addition of a random angle to 
this newly calculated orientation, and we have finally: 
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where is a delta-correlated white noise (^ e [— ^, ^]). 
The interaction introduces a tendency to align (nemati- 
cally) with neighboring particles, so that two simple lim- 
its arise: Complete orientational order settles in the ab- 
sence of noise, whereas particles perform random walks 
for maximal noise (cr = 1). We first characterize the tran- 
sition that necessarily lies in between these two regimes. 
To this aim, we calculate the total order parameter 
Q{N,L) measured for N particles in a square domain 
of linear size L with periodic boundary conditions. We 
use in particular the scalar order parameter 
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which is equal to twice the (positive) eigenvalue of Q, so 
that 5 = 1 for perfect orientational order, and 5* = for 
complete disorder. Starting from random positions and 
orientations, S typically grows in time and eventually 
reaches a statistically stationary state characterized by a 
well-defined distribution function of mean (S). 

Varying the noise intensity a, we observe a continuous 
change of {S). Increasing system size at fixed density 
p — N/L'^, the curves S{N) vs a reveal sharper transi- 
tions for larger sizes, but they do not cross each other 
(Fig.IBi). Figure lit shows that {S{N)) ~ N-i^^l In- 
creasing the noise strength from (7 = towards the tran- 
sition zone, Ci'^) increases from zero to take rather small 
values. Sufficiently deep in the transition zone, the ef- 
fective exponent ({(t,N) = — ^j'^^ can be observed to 
cross over, as ^ oo, from these small values towards 
5 (Fig.^), the value observed at larger a and character- 
istic of a completely disordered phase. All these obser- 
vations are in qualitative agreement with an equilibrium 
KT transition [15|, and signal that only QLRO order is 
present in the ordered phase. At the quantitative level, 
the location of the KT transition point in an equilibrium 
system is characterized hy — jg. Our data is consistent 
with this: for p — ^, we find ^ ~ at CTc — 0.113(5), 
and for larger noise values C(o', N) show signs of crossing 
over to i (i.e. one is in the disordered phase). 

These observations are strengthened by the study of 
the orientational spatial correlation functions g27n{r) = 
(cos[2to(0(O) — 0{r))\) (here m is integer and averages are 
taken both in space and time) . They decay algebraically 
at low noise values {g2m{f) ^ r"™ ^ with rj also increas- 
ing with a) and exponentially in the disordered phase 
(with a diverging correlation length ^ as cr approaches 
the critical point) (Fig. |2K,b). At a quantitative level, 
one expects that in equilibrium 77(17) — 4C((7) in the or- 
dered phase 0]. This is roughly borne out of our data, 
even though good estimates of the correlation functions 
are difficult to obtain close to the transition. In the dis- 
ordered phase, in particular, the expected divergence of 
5 is observed (Fig. Et) but not in its expected functional 
dependence as ^ can only be safely estimated rather far 
away from threshold. Despite these difficulties, we can 
check that the critical noise level determined above is 
consistent with our correlation function data (Fig.[2li). 

The above results are characteristic of QLRO order 
and of a KT-like phase transition. Further preliminary 
results show that the disclination unbinding mech- 
anism characteristic of the KT transition is also at 
work in our non-equilibrium context. All this indicates 
the proper but costly methodology to locate the critical 
noise level Cc- We have not, so far, used this protocol ex- 
tensively, but ongoing simulations indicate that the crit- 
ical line may scale like CTc ~ and extend to arbitrarily 
small densities and noise levels. 

In spite of the equilibrium-like properties of the transi- 
tion to QLRO, the non-equilibrium character of the prob- 
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FIG. 1; (Color online) Transition to nematic order at density 
p = N/L^ = i. (a): (time -averaged) scalar order parameter 
(S) vs noise strength o for various sizes, (b): (S) vs A'^ at 
a = 0.02, 0.04, 0.06, 0.08, 0.1 from top to bottom; the lines 
are fitted power-laws. Each point represents an average over 
10" — 10^ timesteps after transients, (c): same as (b) but 
for larger noise values: a = 0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 
0.17, 0.18. For large a values these curves crossover to 
decay, (d): exponents extracted from the power-laws shown 
in (b) and (c) vs a. The line linking the symbols is a (cubic) 
fit. For the last 2 points, the exponent is only effective, as it 
was estimated from small A'' values, before signs of crossover 
appear. The dotted lines indicate the threshold ctc = 0.113(5) 
estimated from the condition C(''"c) = (see text). 

lem is manifests itself in strong density fluctuations, as 
predicted by Ramaswamy et al. |^ : we measured, in the 
ordered phase, the density fluctuations in square boxes of 
linear size i embedded in a square domain of linear size L. 
These boxes contain, on average, (n) = pi"^ particles. As 
long a,s £ < L, An, the rms of the fluctuations of n, scales 
linearly with (n) (and not \/n), in agreement with 
(Fig.|3J). These giant number fluctuations are the statisti- 
cal consequence of the complex, coupled, spatiotemporal 
dynamics of density and orientation in the system. After 
transients, low- and high-density regions emerge with the 
highly populated domains taking the form of bands in- 
side which nematic order is strong (Fig.^J. These bands 
evolve (move, split, merge, dissolve and form again along 
a new direction) over very long timescales. Typically, 
however, a single band is present at any given time, inde- 
pendently of the system size, and its characteristic evolu- 
tion time grows with system size. Thus, contrary to what 
has been observed in the polar SPP case , these struc- 
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FIG. 2: (Color online) Orientational correlation function 
g2{r) as calculated from the orientation field coarse-grained 
over boxes of linear size 4 (p = ^, time-average over 10^-10^ 
timesteps after transients for each run, L = 256 in (a,d) and 
L = 128 in (b)). (a) in log- log scales at small noise values 
a = 0.02, 0.04, 0.06, 0.08, 0.10, 0.11, 0.12, and 0.13 from top 
to bottom, (b) in lin-log scales for cr = 0.14 to 0.3 from top 
to bottom, (c) variation of the correlation length ^ extracted 
from the exponential tails in (b) as the transition point is ap- 
proached (dashed line at the estimated threshold tJc = 0.113). 
(d) same as in (a), data multiplied by the expected exponent 
at threshold ??c = j: around the estimated threshold value 
(Jc = 0.113(5), the curves are flat and straight. 



tures have no well-defined length- or time-scales. Ongo- 
ing work aims at quantifying these statements |1(tI |. We 
now comment on our results. In our minimal model, 
the density and orientation fields are coupled intimately. 
Whether the mechanism put forward in [6j — namely that 
the current in the conservation equation for the concen- 
tration has a contribution proportional to {dyO, d^O) — is 
actually present here remains to be seen explicitly, but 
we believe that this is the case because the general sym- 
metry arguments invoked there must apply. As a mat- 
ter of fact, in our model, intrinsically non-equilibrium 
features are still observed when the displacement of the 
particles is not made along their orientation — the "nat- 
ural" case if one has in mind particles with an elongated 
physical shape — , but, for instance, perpendicularly to 
their axis (not shown). On the other hand, displacing 
them randomly along one of the four directions defined 
by adding multiples of |- to their current angle yields 
normal fluctuations (An ~ ^/n) and no segregation (not 
shown). This is not surprising since then particles effec- 
tively perform random walks on scales larger than the 
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FIG. 3: (Color online) Root mean square density fluctuations 
An in square boxes containing (n) particles on average (linear 
size ( = ^frifp^ p — ^, time-average over lO^-lO'^ timesteps 
after transients). The saturation at large (n) occurs when 
£ ~ L. Top two curves: giant fluctuations in the ordered phase 
at 2 different system sizes. Bottom curve: normal fluctuations 
in the disordered phase. 



elementary displacement. This case is thus equivalent to 
strictly decoupling density fluctuations from the orien- 
tation fleld by letting particles be non-interacting ran- 
dom walkers. Apart from such "equilibrium" cases, it 
seems that any coupling generically triggers giant fluctu- 
ations, so that they appear as a robust feature, as im- 
licitly implied by the general arguments developed in 
Whereas giant number fluctuations are easily ob- 
served, we have not been able to measure the slow decay 
of t^ged-particle velocity autocorrelations also predicted 
in |y] . We believe this is because the only signiflcant mo- 
tion in our model, apart from the microscopic random 
displacements, is due to the very slow dynamics of the 
high-density bands (Fig. 0)). Thus, we expect such an 
effect to be only observable on timescales so large that 
they are not easily accessible. 

The giant number fluctuations taken as the signature 
of the "non-equilibriumness" of the system are tanta- 
mount to the formation of the high-density ordered band 
described above (Fig. . Whether this is "true" macro- 
scopic phase separation is thus a key question, which 
we just learnt to be the subject of Here we have 

shown that despite this spectacular phenomenon, the 
phase transition is similar, as far as we can tell numeri- 
cally, to the equilibrium one, with an ordered phase char- 
acterized by QLRO only. This has to be paralleled to the 
case of polar SPP, where true long-range order is ascer- 
tained 4], and the transition is discontinuous (i.e. first- 
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FIG. 4: (Color online) Typical time series of S in the ordered 
phase, (p = i, L = 128, a = 0.1). Note the large excur- 
sions to small 5* values which make time-averaging difficult, 
(b-d) snapshots taken during the run shown in (a) at the 
times indicated there. For each particle, a small segment (of 
arbitrary length) centered on its position and aligned on its 
orientation is drawn, (b) and (d) are typical ordered states, 
while (c) represents the more disordered episodes when the 
macroscopic structure changes orientation. 

order like) At equilibrium, the Mermin- Wagner the- 
orem states that true LRO cannot arise from the sponta- 
neous breaking of a continuous symmetry in two dimen- 
sions (18i] . Out of equilibrium, this constraint disappears 
P, Q , and it is quite interesting to notice that polar and 
nematic SPP behave differently. Note also that for polar 
particles the ordered phase is not density-homogeneous, 
but typically consists of well-defined solitary bands with 
high density and strong order moving in a low-density 
disordered background (8||. These bands, which can ap- 
pear in arbitrary numbers depending on the system ge- 
ometry, are very different from the single, fluctuating, 
splitting and merging object described here in Fig. 0] 
Thus, combining the present conclusions and the results 
obtained recently on polar SPP, we see clearly how the 
non-equilibrium nature of these nonlinear driven systems 
can emerge differently at the collective level. 
To summarize, we have introduced a minimal model for 
active nematic particles in which the density and orien- 
tation fields are coupled in a natural way. This non- 
equilibrium model exhibits a KT-like phase transition to 
a nematically-ordered phase in which giant density fluc- 
tuations arise. This constitutes a first non-trivial con- 
firmation of the intrinsically non-equilibrium properties 



predicted in [al and calls for further experimental stud- 
ies, either with assemblies of granular elongated particles 
or, better, in biological systems where such effects could 
play an important role. In view of the notorious difficul- 
ties encountered to decide about similar issues in the case 
of polar particles, our results will need to be confirmed 
by a proper renormalisation group analysis. On the mod- 
eling side, finally, our approach can easily be extended to 
other space dimensions and/or to more complex types 
of interactions. Future work will explore these issues, in 
particular three-dimensional systems and the nature of 
tetratic order out of equilibrium. 

We thank Sriram Ramaswamy for fruitful exchanges 
and the communication of his recent results Il4l . 
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